1 Nood aan een goeie controle

  • Een goede controle groep is cruciaal.

  • Om het effect van een interventie te evalueren is het noodzakelijk om de behandelingsgroep te vergelijken met een controle groep.

  • We moeten ons bewust zijn van Confounding

  • Groepen in een observationele studie zijn dikwijls niet vergelijkbaar. Geavanceerde statistische methodel zijn nodig om causale verbanden te kunnen trekken.

  • Double blinding

  • Gerandomiseerde studies: random toewijzing van subjecten tot de verschillende behandelingsarmen van een studie \(\rightarrow\) vergelijkbare groepen.


2 Randomisatie

  • Randomisatie compleet at random (no systematic allocation).

2.1 Simpele Randomisatie

  • Kan leiden tot verschillen in aantal subjecten per behandelingsarm

  • in 5% of van de steekproeven zien we een imbalans van
    • minstens 60:40 in een studie met 100 subjecten, en
    • van minstens 531:469 in een studie met 1000 subjecten.
  • The imbalans is niet problematisch, maar zorgt voor een verlies in precisie op de schatters (vb. gemiddelden).


2.2 Gebalanceerde Randomisatie

  • Gelijke aantallen van elke behandeling worden toegekend aan een blok van 2, 4 of meer subjecten.
      1. AB, (2) BA
      1. AABB, (2) ABAB, (3) ABBA, (4) BABA, (5) BAAB, (6) BBAA
    • …
  • Gebalanceerde Randomisatie verzekert \(\pm\) gelijk aantal subjecten in elke behandelingsarm van het experiment (vb controle en behandelingsgroep).

  • Zorgt er niet voor dat we een gelijk aantal mannen en vrouwen hebben in elke groep, etc.

  • In kleine studies kunnen de groepen ongebalanceerd zijn in andere karakteristieken (e.g. gender, afkomst, leeftijd …)

  • Dat is niet problematisch omdat het at random gebeurt, maar het zorgt opnieuw voor een verlies in precisie.


2.3 Gestratificeerde Randomisatie

  • De imbalans voor bijvoorbeeld gender kan worden vermeden door gestratificeerde randomisatie: gebalanceerde randomisatie per stratum
Stratified Randomisatie

Stratified Randomisatie


3 Blok Designs

3.1 Gene expression example

  • dm: diabetic medium, nd: non diabetic medium, co: control
  • 4 bio-repeats, 2 technische repeats/biorep

  • dm: diabetic medium, nd: non diabetic medium, co: control
  • 4 bio-reps, 2 techreps/biorep, 2 plates A & B
  • Treatment and plate almost entirely confounded

3.2 Nature methods: Points of significance - Blocking


3.3 Randomised complete block designs

  • Experimentele designs waarbij experimentele subjecten eerst worden opgedeeld in blokken
  • Vervolgens wordt elk niveau van de interventie binnen elk blok toepast en via randomisatie toegewezen.
  • Dat kan b.v. d.m.v. gestratificeerde randomisatie: stratificatie volgens blokken
  • Nuttig wanneer uitkomst sterk varieert tussen subjecten zodat interventie-effect moeilijk op te pikken is
  • Als gegevens veel minder variabel zijn per blok, laat randomiseren van de interventie per blok toe om interventie-effect binnen blok te evalueren met veel minder ruis.
  • Ook veel gebruikt bij technische variabiliteit.
    • Vb grotere experimenten: meerdere celculturen nodig, niet alles kan tegelijk in lab, meerdere sequeneringsruns om genexpressie te meten
    • Fluctuaties in de labo-condities , tussen celculturen of van sequeneringsrun tot sequeneringsrun zorgen voor extra technische ruis.
    • Blokken (vb. tijdstippen, runs, celculturen) en behandelingen Randomiseren binnen elk blok
    • Opnieuw interventie-effecten schatten met veel minder ruis.

4 Steekproefgrootte

  • De steekproefgrootte is ook cruciaal

  • Eenzelfde proefopzet waarbij een steekproef uit dezelfde populatie wordt getrokken maar met meer subjecten \(\rightarrow\) Meer precieze resultaten (Zie introductie)


LS0tCnRpdGxlOiAiMy4gRW5rZWxlIGNvbmNlcHRlbiB2YW4gZXhwZXJpbWVudGVlbCBkZXNpZ24iCmF1dGhvcjogIkxpZXZlbiBDbGVtZW50IgpkYXRlOiAic3RhdE9taWNzLCBHaGVudCBVbml2ZXJzaXR5IChodHRwczovL3N0YXRvbWljcy5naXRodWIuaW8pIgpvdXRwdXQ6CiAgICBodG1sX2RvY3VtZW50OgogICAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCiAgICAgIHRoZW1lOiBjb3NtbwogICAgICB0b2M6IHRydWUKICAgICAgdG9jX2Zsb2F0OiB0cnVlCiAgICAgIGhpZ2hsaWdodDogdGFuZ28KICAgICAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlCi0tLQoKYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9CmtuaXRyOjpvcHRzX2NodW5rJHNldChpbmNsdWRlID0gVFJVRSwgY29tbWVudCA9IE5BLCBlY2hvID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgIG1lc3NhZ2UgPSBGQUxTRSwgd2FybmluZyA9IEZBTFNFKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShOSEFORVMpCmBgYAoKYGBge3IgcG9wMlNhbXAyUG9wLCBvdXQud2lkdGg9JzgwJScsZmlnLmFzcD0uOCwgZmlnLmFsaWduPSdjZW50ZXInLGVjaG89RkFMU0V9CmlmICgicGkiJWluJWxzKCkpIHJtKCJwaSIpCmtvcHZvZXRlcjwtZnVuY3Rpb24oeCx5LGFuZ2xlPTAsbD0uMixjZXguZG90PS41LHBjaD0xOSxjb2w9ImJsYWNrIikKewphbmdsZT1hbmdsZS8xODAqcGkKcG9pbnRzKHgseSxjZXg9Y2V4LmRvdCxwY2g9cGNoLGNvbD1jb2wpCmxpbmVzKGMoeCx4K2wqY29zKC1waS8yK2FuZ2xlKSksYyh5LHkrbCpzaW4oLXBpLzIrYW5nbGUpKSxjb2w9Y29sKQpsaW5lcyhjKHgrbC8yKmNvcygtcGkvMithbmdsZSkseCtsLzIqY29zKC1waS8yK2FuZ2xlKStsLzQqY29zKGFuZ2xlKSksYyh5K2wvMipzaW4oLXBpLzIrYW5nbGUpLHkrbC8yKnNpbigtcGkvMithbmdsZSkrbC80KnNpbihhbmdsZSkpLGNvbD1jb2wpCmxpbmVzKGMoeCtsLzIqY29zKC1waS8yK2FuZ2xlKSx4K2wvMipjb3MoLXBpLzIrYW5nbGUpK2wvNCpjb3MocGkrYW5nbGUpKSxjKHkrbC8yKnNpbigtcGkvMithbmdsZSkseStsLzIqc2luKC1waS8yK2FuZ2xlKStsLzQqc2luKHBpK2FuZ2xlKSksY29sPWNvbCkKbGluZXMoYyh4K2wqY29zKC1waS8yK2FuZ2xlKSx4K2wqY29zKC1waS8yK2FuZ2xlKStsLzIqY29zKC1waS8yK3BpLzQrYW5nbGUpKSxjKHkrbCpzaW4oLXBpLzIrYW5nbGUpLHkrbCpzaW4oLXBpLzIrYW5nbGUpK2wvMipzaW4oLXBpLzIrcGkvNCthbmdsZSkpLGNvbD1jb2wpCmxpbmVzKGMoeCtsKmNvcygtcGkvMithbmdsZSkseCtsKmNvcygtcGkvMithbmdsZSkrbC8yKmNvcygtcGkvMi1waS80K2FuZ2xlKSksYyh5K2wqc2luKC1waS8yK2FuZ2xlKSx5K2wqc2luKC1waS8yK2FuZ2xlKStsLzIqc2luKC1waS8yLXBpLzQrYW5nbGUpKSxjb2w9Y29sKQp9CgpwYXIobWFyPWMoMCwwLDAsMCksbWFpPWMoMCwwLDAsMCkpCnBsb3QoMCwwLHhsYWI9IiIseWxhYj0iIix4bGltPWMoMCwxMCkseWxpbT1jKDAsMTApLGNvbD0wLHhheHQ9Im5vbmUiLHlheHQ9Im5vbmUiLGF4ZXM9RkFMU0UpCnJlY3QoMCw2LDEwLDEwLGJvcmRlcj0icmVkIixsd2Q9MikKdGV4dCguNSw4LCJwb3B1bGF0aW9uIixzcnQ9OTAsY29sPSJyZWQiLGNleD0yKQpzeW1ib2xzICgzLCA4LCBjaXJjbGVzPTEuNSwgY29sPSJyZWQiLGFkZD1UUlVFLGZnPSJyZWQiLGluY2hlcz1GQUxTRSxsd2Q9MikKc2V0LnNlZWQoMzMwKQpncmlkPXNlcSgwLDEuMywuMDEpCgpmb3IgKGkgaW4gMTo1MCkKewoJYW5nbGUxPXJ1bmlmKG49MSxtaW49MCxtYXg9MzYwKQoJYW5nbGUyPXJ1bmlmKG49MSxtaW49MCxtYXg9MzYwKQoJcmFkaXVzPXNhbXBsZShncmlkLHByb2I9Z3JpZF4yKnBpL3N1bShncmlkXjIqcGkpLHNpemU9MSkKCWtvcHZvZXRlcigzK3JhZGl1cypjb3MoYW5nbGUxLzE4MCpwaSksOCtyYWRpdXMqc2luKGFuZ2xlMS8xODAqcGkpLGFuZ2xlPWFuZ2xlMikKfQp0ZXh0KDcuNSw4LCJMZW5ndGUgaW4gcG9wdWxhdGllIixjb2w9InJlZCIsY2V4PTEuMikKCnJlY3QoMCwwLDEwLDQsYm9yZGVyPSJibHVlIixsd2Q9MikKdGV4dCguNSwyLCJzYW1wbGUiLHNydD05MCxjb2w9ImJsdWUiLGNleD0yKQpzeW1ib2xzICgzLCAyLCBjaXJjbGVzPTEuNSwgY29sPSJyZWQiLGFkZD1UUlVFLGZnPSJibHVlIixpbmNoZXM9RkFMU0UsbHdkPTIpCmZvciAoaSBpbiAwOjIpCglmb3IgKGogaW4gMDo0KQp7CgoJa29wdm9ldGVyKDIuMStqKigzLjktMi4xKS80LDEuMStpKQp9CnRleHQoNy41LDIsIkxlbmd0ZSBpbiBzYW1wbGUiLGNvbD0iYmx1ZSIsY2V4PTEuMikKCmFycm93cygzLDUuOSwzLDQuMSxjb2w9ImJsYWNrIixsd2Q9MykKYXJyb3dzKDcsNC4xLDcsNS45LGNvbD0iYmxhY2siLGx3ZD0zKQp0ZXh0KDEuNSw1LCJFWFAuIERFU0lHTiAoMSkiLGNvbD0iYmxhY2siLGNleD0xLjIpCnRleHQoOC41LDUsIlNDSEFUVElORyAmXG5JTkZFUkVOVElFICgzKSIsY29sPSJibGFjayIsY2V4PTEuMikKdGV4dCg3LjUsLjUsIkRBVEEgRVhQTE9SQVRJRSAmXG5CRVNDSFJJSlZFTkRFIFNUQVRJU1RJRUsgKDIpIixjb2w9ImJsYWNrIixjZXg9MS4yKQpgYGAKCi0tLQoKI05vb2QgYWFuIGVlbiBnb2VpZSBjb250cm9sZQoKLSBFZW4gZ29lZGUgY29udHJvbGUgZ3JvZXAgaXMgY3J1Y2lhYWwuCgotIE9tIGhldCBlZmZlY3QgdmFuIGVlbiBpbnRlcnZlbnRpZSB0ZSBldmFsdWVyZW4gaXMgaGV0IG5vb2R6YWtlbGlqayBvbSBkZSBiZWhhbmRlbGluZ3Nncm9lcCB0ZSB2ZXJnZWxpamtlbiBtZXQgZWVuIGNvbnRyb2xlIGdyb2VwLgoKLSBXZSBtb2V0ZW4gb25zIGJld3VzdCB6aWpuIHZhbiBDb25mb3VuZGluZwoKLSBHcm9lcGVuIGluIGVlbiBvYnNlcnZhdGlvbmVsZSBzdHVkaWUgemlqbiBkaWt3aWpscyBuaWV0IHZlcmdlbGlqa2JhYXIuIEdlYXZhbmNlZXJkZSBzdGF0aXN0aXNjaGUgbWV0aG9kZWwgemlqbiBub2RpZyBvbSBjYXVzYWxlIHZlcmJhbmRlbiB0ZSBrdW5uZW4gdHJla2tlbi4KCi0gRG91YmxlIGJsaW5kaW5nCgotIEdlcmFuZG9taXNlZXJkZSBzdHVkaWVzOiByYW5kb20gdG9ld2lqemluZyB2YW4gc3ViamVjdGVuIHRvdCBkZSB2ZXJzY2hpbGxlbmRlIGJlaGFuZGVsaW5nc2FybWVuIHZhbiBlZW4gc3R1ZGllICRccmlnaHRhcnJvdyQgdmVyZ2VsaWprYmFyZSBncm9lcGVuLgoKLS0tCgojIFJhbmRvbWlzYXRpZQoKLSBSYW5kb21pc2F0aWUgY29tcGxlZXQgYXQgcmFuZG9tIChubyBzeXN0ZW1hdGljIGFsbG9jYXRpb24pLgoKIyMgU2ltcGVsZSBSYW5kb21pc2F0aWUKCi0gS2FuIGxlaWRlbiB0b3QgdmVyc2NoaWxsZW4gaW4gYWFudGFsIHN1YmplY3RlbiBwZXIgIGJlaGFuZGVsaW5nc2FybQoKLSBpbiA1JSBvZiB2YW4gZGUgc3RlZWtwcm9ldmVuIHppZW4gd2UgZWVuIGltYmFsYW5zIHZhbgogICAgLSBtaW5zdGVucyA2MDo0MCBpbiBlZW4gc3R1ZGllIG1ldCAxMDAgc3ViamVjdGVuLCBlbgogICAgLSB2YW4gbWluc3RlbnMgNTMxOjQ2OSBpbiBlZW4gc3R1ZGllIG1ldCAxMDAwIHN1YmplY3Rlbi4KCi0gVGhlIGltYmFsYW5zIGlzIG5pZXQgcHJvYmxlbWF0aXNjaCwgbWFhciB6b3JndCB2b29yIGVlbiB2ZXJsaWVzIGluIHByZWNpc2llIG9wIGRlIHNjaGF0dGVycyAodmIuIGdlbWlkZGVsZGVuKS4KCi0tLQoKIyMgR2ViYWxhbmNlZXJkZSBSYW5kb21pc2F0aWUKCi0gR2VsaWprZSBhYW50YWxsZW4gdmFuIGVsa2UgYmVoYW5kZWxpbmcgd29yZGVuIHRvZWdla2VuZCBhYW4gZWVuIGJsb2sgdmFuIDIsIDQgb2YgbWVlciBzdWJqZWN0ZW4uCiAgICAtICgxKSBBQiwgKDIpIEJBCiAgICAtICgxKSBBQUJCLCAoMikgQUJBQiwgKDMpIEFCQkEsICg0KSBCQUJBLCAoNSkgQkFBQiwgKDYpIEJCQUEKICAgIC0gLi4uCgotIEdlYmFsYW5jZWVyZGUgUmFuZG9taXNhdGllIHZlcnpla2VydCAkXHBtJCBnZWxpamsgYWFudGFsIHN1YmplY3RlbiBpbiBlbGtlIGJlaGFuZGVsaW5nc2FybSB2YW4gaGV0IGV4cGVyaW1lbnQgKHZiIGNvbnRyb2xlIGVuIGJlaGFuZGVsaW5nc2dyb2VwKS4KCi0gWm9yZ3QgZXIgbmlldCB2b29yIGRhdCB3ZSBlZW4gZ2VsaWprIGFhbnRhbCBtYW5uZW4gZW4gdnJvdXdlbiBoZWJiZW4gaW4gZWxrZSBncm9lcCwgZXRjLiAgCgotIEluIGtsZWluZSBzdHVkaWVzIGt1bm5lbiBkZSBncm9lcGVuIG9uZ2ViYWxhbmNlZXJkIHppam4gaW4gYW5kZXJlIGthcmFrdGVyaXN0aWVrZW4gKGUuZy4gZ2VuZGVyLCBhZmtvbXN0LCBsZWVmdGlqZCAuLi4pCgotIERhdCBpcyBuaWV0IHByb2JsZW1hdGlzY2ggb21kYXQgaGV0IGF0IHJhbmRvbSBnZWJldXJ0LCBtYWFyIGhldCB6b3JndCBvcG5pZXV3IHZvb3IgZWVuIHZlcmxpZXMgaW4gcHJlY2lzaWUuCgotLS0KCiMjICoqR2VzdHJhdGlmaWNlZXJkZSBSYW5kb21pc2F0aWUqKgoKLSBEZSBpbWJhbGFucyB2b29yIGJpanZvb3JiZWVsZCBnZW5kZXIga2FuIHdvcmRlbiB2ZXJtZWRlbiBkb29yIGdlc3RyYXRpZmljZWVyZGUgcmFuZG9taXNhdGllOiBnZWJhbGFuY2VlcmRlIHJhbmRvbWlzYXRpZSBwZXIgc3RyYXR1bQoKIVtTdHJhdGlmaWVkIFJhbmRvbWlzYXRpZV0oaHR0cHM6Ly9naXRodWIuY29tL3N0YXRPbWljcy9zYmMyMC9yYXcvbWFzdGVyL3RoZW9yeS9maWd1cmVzL3N0cmF0aWZpY2F0aW9uLnBuZyl7IHdpZHRoPTUwJSB9CgotLS0KCiMgQmxvayBEZXNpZ25zCgojI0dlbmUgZXhwcmVzc2lvbiBleGFtcGxlCgotIGRtOiBkaWFiZXRpYyBtZWRpdW0sIG5kOiBub24gZGlhYmV0aWMgbWVkaXVtLCBjbzogY29udHJvbAotIDQgYmlvLXJlcGVhdHMsIDIgdGVjaG5pc2NoZSByZXBlYXRzL2Jpb3JlcAohW10oaHR0cHM6Ly9naXRodWIuY29tL3N0YXRPbWljcy9zYmMyMC9yYXcvbWFzdGVyL3RoZW9yeS9maWd1cmVzL3FwY3JCYWREZXNpZ24xLnBuZyl7IHdpZHRoPTEwMCUgfQoKLSBkbTogZGlhYmV0aWMgbWVkaXVtLCBuZDogbm9uIGRpYWJldGljIG1lZGl1bSwgY286IGNvbnRyb2wKLSA0IGJpby1yZXBzLCAyIHRlY2hyZXBzL2Jpb3JlcCwgMiBwbGF0ZXMgQSAmIEIKLSBUcmVhdG1lbnQgYW5kIHBsYXRlIGFsbW9zdCBlbnRpcmVseSBjb25mb3VuZGVkCgohW10oaHR0cHM6Ly9naXRodWIuY29tL3N0YXRPbWljcy9zYmMyMC9yYXcvbWFzdGVyL3RoZW9yeS9maWd1cmVzL3FwY3JCYWREZXNpZ24yLnBuZyl7IHdpZHRoPTEwMCUgfQoKIyMgTmF0dXJlIG1ldGhvZHM6IFBvaW50cyBvZiBzaWduaWZpY2FuY2UgLSBCbG9ja2luZwoKIVtdKGh0dHBzOi8vd3d3Lm5hdHVyZS5jb20vYXJ0aWNsZXMvbm1ldGguMzAwNS5wZGYpeyB3aWR0aD0xMDAlIH0KCgotLS0KCiMjIFJhbmRvbWlzZWQgY29tcGxldGUgYmxvY2sgZGVzaWducwoKLSBFeHBlcmltZW50ZWxlIGRlc2lnbnMgd2FhcmJpaiBleHBlcmltZW50ZWxlIHN1YmplY3RlbiBlZXJzdCB3b3JkZW4gb3BnZWRlZWxkIGluIGJsb2trZW4KLSBWZXJ2b2xnZW5zIHdvcmR0IGVsayBuaXZlYXUgdmFuIGRlIGludGVydmVudGllIGJpbm5lbiBlbGsgYmxvayB0b2VwYXN0IGVuIHZpYSByYW5kb21pc2F0aWUgdG9lZ2V3ZXplbi4KLSBEYXQga2FuIGIudi4gZC5tLnYuIGdlc3RyYXRpZmljZWVyZGUgcmFuZG9taXNhdGllOiBzdHJhdGlmaWNhdGllIHZvbGdlbnMgYmxva2tlbgotIE51dHRpZyB3YW5uZWVyIHVpdGtvbXN0IHN0ZXJrIHZhcmllZXJ0IHR1c3NlbiBzdWJqZWN0ZW4gem9kYXQgaW50ZXJ2ZW50aWUtZWZmZWN0IG1vZWlsaWprIG9wIHRlIHBpa2tlbiBpcwotIEFscyBnZWdldmVucyB2ZWVsIG1pbmRlciB2YXJpYWJlbCB6aWpuIHBlciBibG9rLCBsYWF0IHJhbmRvbWlzZXJlbiB2YW4gZGUgaW50ZXJ2ZW50aWUgcGVyIGJsb2sgdG9lIG9tICBpbnRlcnZlbnRpZS1lZmZlY3QgYmlubmVuIGJsb2sgdGUgZXZhbHVlcmVuIG1ldCB2ZWVsIG1pbmRlciBydWlzLgotIE9vayB2ZWVsIGdlYnJ1aWt0IGJpaiB0ZWNobmlzY2hlIHZhcmlhYmlsaXRlaXQuCiAgICAtIFZiIGdyb3RlcmUgZXhwZXJpbWVudGVuOiBtZWVyZGVyZSBjZWxjdWx0dXJlbiBub2RpZywgbmlldCBhbGxlcyBrYW4gdGVnZWxpamsgaW4gbGFiLCBtZWVyZGVyZSBzZXF1ZW5lcmluZ3NydW5zIG9tIGdlbmV4cHJlc3NpZSB0ZSBtZXRlbgogICAgLSBGbHVjdHVhdGllcyBpbiBkZSBsYWJvLWNvbmRpdGllcyAsIHR1c3NlbiBjZWxjdWx0dXJlbiBvZiB2YW4gc2VxdWVuZXJpbmdzcnVuIHRvdCBzZXF1ZW5lcmluZ3NydW4gem9yZ2VuIHZvb3IgZXh0cmEgdGVjaG5pc2NoZSBydWlzLgogICAgLSBCbG9ra2VuICh2Yi4gdGlqZHN0aXBwZW4sIHJ1bnMsIGNlbGN1bHR1cmVuKSBlbiBiZWhhbmRlbGluZ2VuIFJhbmRvbWlzZXJlbiBiaW5uZW4gZWxrIGJsb2sKICAgIC0gT3BuaWV1dyBpbnRlcnZlbnRpZS1lZmZlY3RlbiBzY2hhdHRlbiBtZXQgdmVlbCBtaW5kZXIgcnVpcy4KCi0tLQoKIyBTdGVla3Byb2VmZ3Jvb3R0ZQoKLSBEZSBzdGVla3Byb2VmZ3Jvb3R0ZSBpcyBvb2sgY3J1Y2lhYWwKCi0gRWVuemVsZmRlIHByb2Vmb3B6ZXQgd2FhcmJpaiBlZW4gc3RlZWtwcm9lZiB1aXQgZGV6ZWxmZGUgcG9wdWxhdGllIHdvcmR0IGdldHJva2tlbiBtYWFyIG1ldCBtZWVyIHN1YmplY3RlbiAkXHJpZ2h0YXJyb3ckIE1lZXIgcHJlY2llemUgcmVzdWx0YXRlbiAoWmllIGludHJvZHVjdGllKQoKCi0tLQoKIyBbSG9tZV0oaHR0cHM6Ly9ndHBiLmdpdGh1Yi5pby9QU0xTMjAvKSB7LX0K